Temperature and magnetic-field dependence of the conductivity of YBa 2 Cu 3 7 _<j films 
in the vicinity of superconducting transition: Effect of Tc-inhomogeneity 
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Temperature and magnetic field dependences of the conductivity of YBa2Cu307_f films in the 
transition region are analyzed taking into account spatial inhomogeneity in transition temperature, 
T c . (i) An expression for the superconducting contribution to conductivity, <j s (T, H, T c ), of a homo- 
geneous superconductor for low magnetic fields, H <C H C 2(0), is obtained using the solution of the 
Ginzburg-Landau equation in form of perturbation expansions [S.Ullah, A.T.Dorsey, Phys. Rev. B 
44, 262 (1991)]. (ii) The error in a s (T, H,T C ) occurring due to the presence of T c - inhomogeneity 
is calculated and plotted on an H — T plane diagram. These calculations use an effective medium 
approximation and a Gaussian distribution of T c . (iii) Measuring the temperature dependences of 
a voltage, induced by a focused electron beam, we determine spatial distributions of the critical 
temperature for YBa2Cu307_« microbridges with a 2 /im resolution. A typical redistribution dis- 
persion is found to be ~1 K. For such dispersion, error in a 3 (T, H,T C ) due to T c -inhomogeneity 
exceeds 30% for magnetic fields H < 1 T and temperatures \T — T c \ < 0.5 K . (iv) Experimental 
R(T, H) dependences of resistance are well described by a numerical solution of a set of Kirchoff 
equations for the resistor network based on the measured spatial distributions of T c and the expres- 
sion for cr s (T, H,T C ). 

PACS numbers: 74.62.-c, 74.25.Fy, 74.80.-g, 74.72.Bk 



I. INTRODUCTION 

The complicated crystal structure of high-T c supercon- 
ductors (HTSC) leads to their substantial spatial inho- 
mogeneity, which is especially important because of the 
very short coherence length, £. Inhomogeneities with 
spatial scale much larger than £ allow for an inhomoge- 
neous distribution of the critical temperature, T c , which 
affects properties of HTSC in the superconducting tran- 
sition vicinity. As a result, it is often difficult to establish 
whether observed behavior of superconductors arise from 
intrinsic properties or from spatial inhomogeneity. This 
impedes analysis of experimental data in the transition 
region, which is often used to determine microscopic su- 
perconducting parameters and the mechanism of super- 
conductivity. 

The most obvious origin of T c -inhomogeneity is 
variation in oxygen content over the sample. For 
YBa2Cu307_,5 (YBCO), T c is a relatively weak function 
of S at 6.85 < 7 - S < 7 Lsa called 90 K plateau) and 
falls abruptly at higher Sun Even S variations within 
the 90 K plateau can lead to w 1 K variations in 
T c . Meanwhile, experimental x-ray data show that 5 
variation can be substantially higher even for crystals 
exhibiting excellent transport properties a Another ori- 
gin of T c -inhomogeneity is variation in cation (Y, Ba, 
Cu) composition. This origin can be dominant in thin 
YBCO films, as shown by simultaneous spatially re- 
solved studies of cation composition and T c using elec- 
tron probe microanalysis and low-temperatur.e|-£cannmg 
electron microscopy (LTSEM), respectively.crcl Elastic 



stresses around structural defects can also lead to Tc- 
inhomogeneity due to a .-strong pressure dependence of 
T c in HTSC compounds.Q T c increase due to edge dislo- 
cations awi low-angle grain boundaries was calculated to 
be ~ 1 K.I3 Spatial variations of the c-axis lattice parame- 
ter revealed by x-ray studies in YBCO films with almost 
uniform oxygep. content also suggest stress-induced T c - 
inhomogeneityjj 

Presence of Tc-inhomogeneity manifests itself in vari- 
ous HTSC properties. Temperature dependence of the 
depinning current density in YBCO crystals implies, that 
pinning sites are induced by spatial variations of T c Ja Sys- 
tematic studies of YBCO crystal magnetization curves 
suggest the presence of local regions with reduced rOKpi 
gen content, leading to the so called peak effect .E2rO 
Meanwhile, Tc-inhomogeneity should have even greater 
impact on temperature dependences of transport coeffi- 
cients just above the superconducting transition. This is 
confirmed by experimental data on conductivity, magne- 
toconductivity, and the Hall coefficient at temperatures 
T > T c + 2 K from Refs. [l3]-|!I] wm ch were explained 
by assuming a Gaussian distribution of T c with disper- 
sion in the range 0.6-2.3 K. However, these results can 
only serve as an indirect indication of the presence of Tc- 
inhomogeneity, due to a lack of experimental data about 
real T c -distribution in the samples. Moreover, the tem- 
perature region in the vicinity of the superconducting 
transition, where Tc-inhomogeneity is especially impor- 
tant, was not considered. 

A step forward has been made in Ref. [l^ where resis- 
tor network calculations are used to analyse current den- 
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sity redistributions in T c -inhomogeneous supercondictor 
in the transition region. It is shown that some anomal- 
ities in the temperature dependence of in-plane magne- 
toresistivity, such as negative magnetoresistivity excess, 
which are usually attributed to intrinsic effects can be 
quantitatively expalined by nonuniform T c -distribution. 

In the present work, we investigate the influence of 
T c -inhomogeneity on properties of HTSC throughout the 
transition region and analyze experimentally determined 
spatial distributions of T c (T c -maps). Measuring T c -maps 
of YBCO films by LTSEM with 2 /im resolution, we re- 
veal k 1 K scatter of T c over the films. To calculate the 
effective conductivity of such an inhomogeneous material, 
one needs the expression for conductivity <r(T, H, T c ) of a 
uniform superconductor valid throughout the transition 
region. Such an expression was obtained in Ref. [l7| by 
solving the time-dependent Ginzburg-Landau equation 
with Lawrence and Doniach Hamiltonian in the Hartree 
approximation. 

It is well-known that magnetic field leads to a broad- 
ening of the superconducting transition which, is roughly 
proportional to (dH c2 /dT)~l To « 0.5 K/TEJ For fields 
H > 2 T, this broadening dominates over inhomogeneous 
broadening due to scatter of local values of T c . There- 
fore, we are interested in low fields, H < 2 T, where the 
influence of T c -inhomogeneity is essential. Moreover, this 
range of magnetic fields is actual for most HTSC applica- 
tions. Unfortunately, the final formula for conductivity 
obtained in Ref. |l7| is only valid for high magnetic fields. 
In the present work, we deduce expressions for conduc- 
tivity valid for low magnetic fields, H -C H C 2, from the 
perturbation expansions in Ref. [Tt]. 

The present paper is organized as follows. In Sec. II the 
expression for the Cooper pair conductivity, <r s (T, H, T c ), 
of a homogeneous superconductor is derived. In Sec. Ill 
we discuss methods to calculate the effective conductiv- 
ity of an inhomogeneous superconductor. In Sec. IV the 
error in the value of a s (T, H,T C ) occurring due to the 
presence of T c -inhomogeneity is calculated and the re- 
sults are plotted on the H — T plane diagram. Sec. V 
describes the samples and experimental techniques. In 
Sec. VI, measured spatial distributions of T c are ana- 
lyzed and the broadening of the superconducting transi- 
tion due to Tc-inhomogeneity is discussed. Finally, the 
experimental R{T, H) dependences of resistance are in- 
terpreted on the basis of measured T c -distributions and 
the expression for cr s (T, H,T C ) derived in Sec. II. 



II. CONDUCTIVITY OF A HOMOGENEOUS 
SUPERCONDUCTOR 

To describe the temperature dependence of conduc- 
tivity of a homogeneous superconductor throughout the 
transition region we employ the results obtained by 
Ullah and Dorsey.tZl They studied the time-dependent 
Ginzburg-Landau equation for anisotropic superconduc- 
tor with the Hamiltonian introduced by Lawrence and 



Doniacht3 and an additional noise term. The magnetic 
field H was assumed to be applied along the c-axis. Us- 
ing the Hartree approximation Ullah and Dorsey ob- 
tained expressions for the transport coefficients which 
gave smooth interpolation between the high-temperature 
regime dominated by Gaussian fluctuations and low- 
temperature flux-flow regime. The expression for the 
Cooper pair conductivity in linear order to electric field 
was obtained in the form of two coupled perturbation 
expansions!!! 

N 

(t s = to Y,( n + + K+i - iA n+1/2 ), (l) 

71=0 

whereE] 

A n = A n (i H , h) = [{e H + 2hn)(l + d 2 (~e H + 2hn))] ~ 1/2 , 

(2) 

and 

N 

e„ = e H -nThJ2 A rv (3) 

n=0 

Here h = H/H C 2(Q), and en is a field-dependent dimcn- 
sionless temperature, en = T/T c — 1 + h. Further, 



n 



8tt 2 (2k 2 - 1) £c(0)fci 
7^ 



(4) 



£,ab(0) and £ c (0) are the correlation lengths in the CuO 
plane and transverse to it: £ o b(0) = %/ y / 2m a i,ao (with 
similar relation for £ c (0)); m a b is the Cooper pair mass 
in CuO plane, ao is related to the parameter a in 
Ginzburg-Landau Hamiltonian as a — ao (T/T c — 1), 
d = s/2£ c (0), where s denotes spacing between CuO 
planes, 7 = Cc(0)/Cab(0) is an anisotropy parameter, (f>o is 
the flux quantum, k is the Ginzburg-Landau parameter, 
H C 2(0) = </)o/27r^ 2 b (0), N — 1/h, and ao is a constant 
with dimensionality of a conductivity. 

In order to avoid summation in the above expressions, 
the following approximation suggested in Ref. 17 is usu- 
ally used. For high magnetic field (in *C 2h) and 3D case 
(d 2 en *C 1) only terms containing Aq are left in Eqs. (|l|) 
and (||) yielding a s = (Jq/^/Ih and e H = e H - Qh/y/cR 
respectively. This leads to 
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where function T{x) satisfies cubic equation, 

xT 2 = \-T' i . (6) 
The solution of this equation can be written down as 
T(x) = 6 + x 2 /(99) -x/3, 

6 = (l/2 - x 3 /27 + v/(27-4.x 3 )/108) V3 . (7) 
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The function T(x) is equivalent to the function F 3 d in 
Rcf. [I?] though expressions (g) and (0) for this function 
are not presented there. 

Despite Eq. (||) is widely used, its applicability range 
needs a special discussion. Indeed, since the series @ are 
diverging, omitting all terms except the first one is hardly 
permissible. At least, it is obviously incorrect if the con- 
dition £h "C 2h does not hold. Hence, for magnetic fields 
we deal with, H < 2 T, the high-field approximation of 
Ref. [H] is not valid except for the low-temperature part 
of the superconducting transition. As argued in Ref. [Tt], 
Eq. H can also be considered as a scaling relation with 
unknown scaling function T and then it is valid in a 
wider range of magnetic fields. However, simple numeri- 
cal calculations show that, e. g., for YBa2Cu307_5, the 
scaling does not work for fields H < 2 T. Therefore, 
we derive new, low-field approximation for conductivity 
from Eqs. (|l|) and (||). For h <C 1 one can replace sum- 
mation .for n > 1 by integration using Euler-Maklaurin 
formula! 2 -!! and obtain the following analytical expressions 
for conductivity: 



a s /a Q =A - 2A 1/2 + -A x - 2A 3/2 + A 2 



3 + 8d 2 3h 
' 4(2 + Ad 2 ) 3 / 2 + "8 



^(l + 2d 2 (e„ + 2h))Al (8) 




(9) 



Equations (^) and (^J) present non-explicit dependence 
of conductivity for a homogeneous superconductor on 
temperature and magnetic field. They are derived with- 
out any assumptions about 3D or 2D character of su- 
perconductivity and therefore applicable for arbitrary 
anisotropy parameter. These equations are used below 
for calculations presented in Sec. Ill and V. 

Under conditions 



e#,ft< 1/d 2 < 1 
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equations and (J9|) can be substantially simplified 
yielding an explicit expression for the Cooper pair con- 
ductivity: 
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7h + 2i H 



where T is the function defined by Eqs. (g) or (|7|). An- 
other advantage of these equations is that the conductiv- 
ity depends on d only through parameter A which charac- 
terizes the shift of the apparent transition temperature 
with respect to T c . This shift is always present in the 
Hartree approximation as a result of renormalization of 
the parameter a in the Ginzburg-Landau Hamiltonian. 



Eqs.(|) and <g) as well as Eqs.Q and @ defining 
temperature and magnetic field dependence of conductiv- 
ity contain only two key parameters: d and Q. Parameter 
d depends on the value of coherence length £ c (0) which 
is not known well. An accurate determination of £ c (0) 
is difficult because of bilayered structure of the YBCO 
unit cell. A systematic analysis performed in Ref. |22| 
shows that a good approximation for YBCO is the as- 
sumption of equally spaced CuO layers with interlayer 
distance s — 6 A. Then, a simultaneous fitting of con- 
ductivity, magnetoconductivity and susceptibility data 
gives £ c (0) = 1.2l leading to d = 2.5B- 
not contradict to results of other worksE 
the distance between superconducting layers s sa 11.7 A. 
They give values of £ c (0) in the range between 1.3 and 
3 A which corresponds to 1.9 < d < 4.5. Obviously, this 
uncertainty in d is very large. Fortunately, it is not so im- 
portant in the range of H and T that we consider. Indeed, 
for YBCO the conditions ([l0]) are satisfied within several 
kelvins around the transition unless the magnetic field is 
very high. Then, as follows from Eqs. ( jll"| ) and (12), 
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i H = {VlTh) 2 ' 3 T- 



(1 + A) T/T c -l + h 



{riTh)V 3 



a = nr r 



ln8ri 2 
2d 



2::.: 
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100 - 




(12) 



FIG. 1. Temperature dependence of Cooper pair conduc- 
tivity for YBCO in magnetic fp)d H = 2 T. Symbols show 
the exact result of UD modelJiJ Eqs. ([!]), and (^); dashed 
line is the high-field approximation proposed by Ullah and 
Dorsey, Eqs. (H) and tXm (shifted along x-axis by —A); solid 
line is the low-field approximation proposed in the present pa- 
per, Eqs. (^) and (^; dotted line is an approximate explicit 
expression for low fields given by Eqs. ([n]) and (|l2]). 
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change in d leads only to a shift of apparent transition 
temperature but does not affect the shape of the tran- 
sition curve. The second parameter f2 depends on the 
coherence length £ a b(0) and the Ginzburg-Landau pa- 
rameter k. The former was also taken from Ref. |2^, 
£ a b = 11 A. Meanwhile, it is quite difficult to find in 
literature an accurate estimate for n. We used the value 
k = 30 providing the best fit to our experimental data 
for YBCO films. Then, from equation (^) one obtains 
!l = 5x 1CT 4 . 

Figure [l] presents comparison of the temperature de- 
pendences of Cooper pair conductivity given by straight- 
forward summation in ([!]) and (^), and by two analytical 
approximations: high-field approximation, Eqs. da) and 
([?]) , proposed by Ullah and Dorsey and low-field approx- 
imation, Eqs. (||) and (||), suggested in the present work. 
An approximate explicit expression for low fields given 
by Eqs. (|ll|) and ( |l2] ) is also shown as dotted line. For 
H = 2 T which is the case shown in the figure, the low- 
field approximation is far more accurate than the high- 
field one. For lower magnetic fields the deviation be- 
tween the result of exact summation and the low-field ap- 
proximation is almost indistinguishable. By contrast, the 
high-field approximation fails for temperatures T T c 
where it four times overestimates the result of exact sum- 
mation which is er = oq /{A^/Eh) . 

Figure clearly shows that the apparent transition tem- 
perature is shifted downward from T c . For given set of 
parameters the dimensionless shift is A w 0.03. In or- 
der to avoid confusion, the data in all figures below are 
shifted along £-axis so that T c corresponds to the appar- 
ent transition temperature. It should be also noted that 
the shift A does not enter the high-field approximation 
suggested by Ullah and Dorsey, Eqs. (||) and This 
approximation predicts transition at T = T c in contra- 
diction to basic equations of UD model, Eqs. ([!]) and ||). 
In order to make the high-field approximation merge all 
other curves in Fig. [j]at least at low T, we had to shift the 
corresponding dashed curve on the value A "by hand" . 

The constant do entering UD model depends on a phe- 
nomenological quantity, the relaxation rate of the order 
parameter. It is natural la estimate oo using well-known 
Aslamazov-Larkin resulted for high-temperature asymp- 
totic in 3D case: a 3 s D = e 2 /32fr£(0)^i". Thus, we have 
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Let us now discuss the applicability range for 
the results obtained in this Section. The indirect 
(Maki-Thora»son) contribution to the order parameter 
fluctuationscJ is not taken into account in the UD model. 
However, there are grounds to believe that neglecting 
Maki-Thompson term would not affect the results ob- 
tained in the transition region (few Kelvins around T c ) 
since the direct Aslamazov-Larkin process is dominant 
over the indirect one in this temperature range.EH One 
should also keep in mind that the UD model does not 
take into account vortex pinning and predicts flux-flow 



behavior in the limit of low temperatures. Therefore, it 
cannot be used at temperatures well below T c , where the 
current- voltage characteristics are nonlinear. 



III. ACCOUNT OF r c -INHOMOGENEITY 

Let us now consider how the properties of a supercon- 
ductor can be affected by spatially inhomogeneous dis- 
tribution of critical temperature. First, we suppose that 
the correlation length r c of T c -distribution is so large that 
the temperature region near T c where £(T) > r c can be 
ignored. This assumption seems to be quite reasonable 
since the coherence length of HTSC is much smaller than 
r c obtained from LTSEM data, see Tab. 1. The condition 
r c 3> £ makes it possible to ignore the correlation be- 
tween the superconducting order parameter in adjacent 
fragments and to consider them independently. There- 
fore, the expression for the conductivity obtained in the 
previous section for a homogeneous superconductor can 
be used to describe local conductivity cx(T, H, T c ) of a 
homogeneous fragment with given T c . 

The straightforward way to determine the conductivity 
a (T, H) of a T c -inhomogeneous superconductor is to 
start from the spatial distribution of T c over the sample. 
The value of tr mh (T, H) can be determined exactly from 
the values of local conductivities tr(T, H,T C ). In this work 
we determined the spatial distributions of the critical 
temperature in YBCO films using LTSEM (see Sec. V). 
This method, however, leads to the lack of information 
about small-scale inhomogeneities with r c < r cxp , where 
r eX p is the spatial resolution of the technique. Therefore, 
if small-scale inhomogeneity is essential, or if spatial dis- 
tribution of T c is unknown, a Gaussian T c -distribution 
function together with, e. g., effective medium approach 
can be used to find cr inh (T, H). 

The problem of conductivity of an inhomogeneous 
medium has the exact analytical solution only for a 
special .case of symmetric distribution of phases in 2D 
system.Ea In the general case one has to use some appro*-, 
imation. According to the effective medium approached 
(EM A), the conductivity cr inh (T, H) is given by the solu- 
tion of the equation 



<J inh (T, H)-a(T, H,T C ) 
(D - 1) cr inh (T, H) + cr(T, H, T c 



-f(T c )dT c = 0, (14) 



where D is the dimensionality of the system. Here f(T c ) 
is a distribution function of critical temperature over the 
sample which shows the relative volume occupied by frag- 
ments with given T c . Despite the apparent simplicity, 
EMA gives rather high accuracy (up to few percents) un- 
less the svstem is in the very vicinity of the percolation 
thresholdB3 In the case of thin film samples with thick- 
ness less than the correlation length of T c -inhomogeneity, 
r c , EMA expression ( |l4| ) with dimensionality D=2 should 
be used. It should be emphasized that this dimension- 
ality has nothing to do with the dimensionality of the 
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superconducting properties mentioned in relation with 
formula (||); the first one depends on the geometry of the 
sample, while the latter is associated with anisotropy of 
the crystal structure. 



IV. H-T DIAGRAMS 

In this section we estimate the effect of T c - 
inhomogeneity on the apparent value of the Cooper 
pair conductivity in the vicinity of the superconduct- 
ing transition. Usually, experimental data on <x(T, in- 
dependences in the transition region are studied first 
by subtracting the conductivity of normal electrons, a n , 
and then analyzing the remaining conductivity of Cooper 
pairs, a s . In the case of inhomogeneous sample such pro- 
cedure would lead to an error in a s : its apparent value 
determined from experimental data would be different 
from that for a homogeneous superconductor. To quanti- 
tatively estimate this error we consider two samples: uni- 
form, with critical temperature T c o, and inhomogeneous 
one with a Gaussian distribution of critical temperatures 
with average T c q and dispersion ST C : 



f(T c ) 



1 



exp 



(r c - T c0 y 

2 5T? 



(15) 



Now two quantities, a\ lom (T,H) and af h (T,H) can be 
compared. <T^ om is the Cooper pair conductivity for ho- 
mogeneous sample given by the expressions (||) and (||) 
obtained on the basis of UD model. The conductiv- 
ity er ml1 of the inhomogeneous sample is determined by 
EMA formula dTJ) with f(T c ) being Gaussian distribu- 
tion function (|l5|), and local conductivities defined as 
sum of the superconducting, <r' lom , and normal, cr„, con- 
tributions. Then, one should subtract the normal contri- 
bution from cr mh and obtain the apparent superconduct- 
ing contribution to conductivity for the inhomogeneous 
sample: 



< lh (T, H) = (j mh (T, H) - a n (T, H). 



(16) 



Further, to proceed with calculations some assump- 
tions are needed about the temperature and magnetic 
field dependences of a n . We neglect the magnetoresis- 
tance of HTSC in the normal state which is very small 
and use a linear approximation for the temperature de- 
pendence of the resistivity: 



a n (T, H) = a n (T, 0) = 1/ (Ci + C 2 T) . 



(17) 
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FIG. 2. Diagrams in the "magnetic field" — "temperature" plane illustrating the effect of T c -inhomogeneity on the apparent 
value of Cooper pair conductivity, a s . Absolute values of the reduced difference of (a) conductivities |oi, nh / 'a^ oul — 1| and (b) 
magnetoconductivities | (a" lh ) ff / ( a ^ 0rn ) H ~~ -M f° r homogeneous and inhomogeneous superconductors are shown. Calcula- 
tions are based on Eqs. Q and (g). For inhomogeneous superconductor, the calculations use an effective medium approach, 
Eq. (^), and a Gaussian distribution of T c with dispersion 8T C = 1 K. Brighter regions correspond to a stronger influence of 
Tc-inhomogeneity. The influence is maximal in low magnetic fields and in the vicinity of T c . Presented results are not valid 
below the white dashed lines which correspond to the melting transition of vortex lattice (the data are taken from Refs. pl|,p2|). 
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The key parameter for calculations is the dispersion of 
Gaussian T c -distribution (H), ST C . We used the value 
ST C = 1 K which is approximately the average dispersion 
for studied YBCO films determined from their T c -maps. 
The values for £1 and d were the same as in Sec. II, other 
parameters were: Ci = 1.06 jtxOm, (72=0.0035 /zfim/K 
(see Tab. 1), T c =90 K, D=2. We also assume that the 
inhomogcncity of a superconductor manifests itself in in- 
homogeneity of the critical temperature only, while all 
the other superconducting parameters and the normal 
state conductivity are supposed to be uniform. 

It is convenient to consider H—T diagram which shows 
the absolute value of the relative difference: {af 11 / a^° m — 
1|, see Fig. ||a. The effect of Tc-inhomogeneity on the 
magnetoconductivity is illustrated by Fig. 0b showing the 
same diagram for the quantity | (<r™ h ) „ / (ffj° m ) H — 1|, 
where (cr s )' H denotes partial derivative of conductivity 
with respect to magnetic field. Brighter regions on the 
diagrams correspond to larger values, i.e., to stronger 
influence of T c -inhomogeneity on the values of <j s and 
(a s )' H . The influence becomes crucial in a 1 K-wide re- 
gion around T c and for magnetic fields H < 0.5 T where 
ignoring T c -inhomogeneity would lead to ~ 50% error in 
a s . The following conclusions can be drawn from the 
diagrams: 

(i) Tc-inhomogeneity plays greater role in the very 
vicinity of the transition; far from the transition the dif- 
ference in local T c 's is small compared to \T — T c \ and, 
hence, not so important. 

(ii) T c -inhomogeneity plays greater role in low mag- 
netic fields. The application of magnetic field leads to 
a broadening of the transition even in a homogeneous 
superconductor. Sipcc for most HTSC dH&jdT sw 
—2 T/K at T = T c ,Ea one can roughly estimate the in- 
crease in the transition width as one degree for H increase 
of 2 T. Therefore, for fields H > 2 T the dispersion in 
critical temperatures 5T C rj 1 K is masked by -ff-induced 
broadening of the transition. 

(iii) Tc-inhomogeneity has greater effect on the magne- 
toconductivity of a superconductor than on its conduc- 
tivity. From practical point of view it is often prefer- 
able to analyze experimental data on magnetoconduc- 
tivity rather than on conductivity. This is because the 
contribution of normal electrons to magnetoconductivity 
is negligible in the vicinity of T c , while the analysis of 
conductivity data always requires account of the normal 
conductivity and, hence, additional assumptions about 
its temperature dependence. However, as follows from 
the diagrams, the analysis of magnetoconductivity data 
needs more careful account of T c -inhomogeneity. Th© 
reason for that, as was earlier noted by Lang et al.fB 
lies in the stronger dependence of magnetoconductivity 
on T c , e.g., for high temperatures, T T c , one has 
o s oc (T — T c )~ 1/2 , while da s /dH oc (T — T c )~ 3/2 . 

The dashed line in Fig. corresponds to the melting 
transition of the Abrikosov vortex l attic e as determined 
from experiments on YBCO crystals.ElIca It is remarkable 



that different methods, neutron small angle scattering,^, 
as well as magnetization and transport measurements,^ 
yield the same position of the melting line. We believe 
that it can serve as a rough estimate of the applicabil- 
ity range of the UD model. Below this line, our results 
obtained on the basis of the UD model are not valid. 



V. SAMPLES AND EXPERIMENTAL DETAILS 

YBa2Cu307_5 films with thickness of 0.2 /im were 
grown by dc magnetron sputtering on NdGaC>3, AlLa03 
and MgO substrata^ The details of the procedure are 
described elsewhere.H X-ray data have shown the pres- 
ence of only (001) reflexes confirming c-orientation of the 
films. The Raman spectroscopy analysis has revealed 
their high epitaxiality. Microbridges of 500x50 fim size 
were formed by a standard photolithography. Six sam- 
ples were investigated; some important parameters are 
presented in Tab. 1. 

The temperature dependences of the resistivity were 
measured at driving current 1 mA and magnetic fields 
H=0, 0.3, 0.6, 0.9 T applied along the c-axis. Mea- 
surements were done inside a temperature stabilized Ox- 
ford He flow cryostat (model CF-1200) under helium at- 
mosphere, using the standard four-probe dc method, a 
Kcithly 220 programmable current source and a Keithly 
182 sensitive digital voltmeter. Contacts to the sam- 
ples were made by thin gold wires attached to the sam- 
ple surface by silver paste. The temperature inside the 
cryostat was controlled and stabilized by an Oxford pro- 
grammable temperature controller ITC4 with accuracy 
up to 0.01 K. The temperature of the sample was mea- 
sured by copper-constantan thermocouple; voltage was 
read by a DMM5000 integrating digital multimeter. The 
measurement was started when the sample was in the 
normal state (at least 40 K above T c ) and performed 
during a slow cooling procedure down to zero resistiv- 
ity of the sample. Then the sample was heated and the 
measurement was repeated at another value of magnetic 
field. The accuracy of the voltage measurements was 
about 10 nV. 

The LTSEM measurements were carried out with an 
automated scanning electron microscope CamScan Series 
4-88 DV100. The microscope is equipped with a cooling 
sample stage, its temperature can be lowered down to 
77 K using an Oxford N flow cryostat. The temperature 
is maintained in the range 77-350 K With accuracy up 
to 0.1 K by a temperature controller ITC4. The bias 
current was varied from 0.2 to 2.0 mA so that its value 
was large enough to detect electron beam induced voltage 
(EBIV) and small enough to avoid distortion of the su- 
perconducting transition. EBIV was measured using the 
standard four-probe method. A precision instrumenta- 
tion amplifier incorporated into the microscope chamber 
was used to increase the signal-to-noise ratio. To extract 
the local EBIV signal, lock-in detection was used with a 
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beam-modulation frequency of 1 kHz. The electron beam 
current was 10 -8 A, while the acceleration voltage was 
10 kV. 

The method for determination of the spatial dis- 
tribution j^aij-critical temperature is based on LTSEM 
techniquetJ'ta and is described in detail in Ref. |37]. Heat- 
ing by electron beam elevates the temperature locally by 
<5Theat < 1 K causing a change, dp, in the local resistiv- 
ity. As a result, a change in the voltage, EBIV, occurs 
across the sample biased by a constant transport current. 
Temperature dependence of EBIV has the maximum at 
some temperature, T m , corresponding to the maximum 
in Sp. Thus, the local transition temperature can be de- 
termined as T c = T m + (5Th C at/2. Scanning the electron 
beam over the film allows us to determine the spatial dis- 
tribution of T c . In order to remove the distorting effects 
associated with thermal diffusion into adjacent regions 
of the film a numerical deconvolution method was used. 
Figure || shows temperature dependences of EBIV for 
two adjacent regions of sample 1 before and after the 
deconvolution procedure. After the deconvolution, both 
dependences have a pronounced major peak; its position 
defines the local T c . It follows from Fig. || that the differ- 
ence in T c for two regions separated by 5 pm can be as 
large as 0.7 K. The method allows the spatial resolution 
of 2 pm and the temperature resolution of 0.2 K. The 
T c -map for sample 1 is shown in Fig. |]. After T c -map is 
determined, one can easily calculate distribution function 
f(T c ) defining the relative volume occupied by fragments 
with given T c ; f(T c ) for sample 1 is shown in Fig. ^| 



Further, using T c -map and the expression for conduc- 
tivity cr(T, H, T c ) of a T c -uniform fragment, one can cal- 
culate the spatial distribution of current density in the 
superconductor. First, the film is approximated by a 
square network of resistors. Then, the set of Kirchoff 
equations is solved with respect to electric potentials in 
the nodes of the network. For this purpose an iterative 
procedure with overrelaxation method is used with-fipted 
potentials of the two opposite sides of the network.E3 As 




FIG. 4. Spatial distribution of critical temperature in 
sample 1 determined from LTSEM data. Distribution is 
smoothed with respect to the initial one measured with 2 pm 
resolution. 
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FIG. 3. Temperature dependences of electron beam in- 
duced voltage measured by LTSEM for two regions of sam- 
ple 1 separated by 5 pm distance. The upper panel shows 
raw signals, the lower panel shows the same dependences after 
deconvolution procedure. Local values of the critical temper- 
ature, Tci and T C 2, are determined by positions of the peaks. 
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FIG. 5. Experimental temperature dependence of resistiv- 
ity (circles) for H = and the distribution function, f(T c ), 
of critical temperature (solid line) determined from T c -map 
shown in Fig. Q All the data are for sample 1. The dotted 
line shows a plausible shape of the total distribution function. 
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a result, the current density distribution as well as the 
total resistance of the superconductor are calculated for 
any temperature in the vicinity of the superconducting 
transition. 

As the temperature is lowered, the current density 
distribution becomes noticeably inhomogeneous. As a 
result, some normal-conducting regions of the film are 
shunted by surrounding superconducting regions. T c of 
these shunted regions cannot be measured by the present 
method. However, one can expect that ambiguity in 
their T c 's would not lead to substantial errors in results 
of resistor network calculations. Indeed, in the high- 
temperature part of the superconducting transition, the 
conductivity of these regions is known since they are in 
the normal state, while at lower temperatures they are 
off the main current path and make a minor contribution 
to the film resistance. Figure || represents distribution 
function f(T c ) determined from the T c -map. The dot- 
ted line in the same figure shows a plausible shape of the 
total distribution function. 

When the spatial distribution of critical temperature 
T c (r) is given, one can estimate the correlation length 
of T c -inhomogeneity. It is defined from the correlation 
function G(r) of ^-distribution : 

= T^R+jT.fR)-^ 

T 2 - T c 

where averaging is performed over all R and all directions 
of r within the bridge. The value G = 1 corresponds to 
full correlation and G = to the absence of correlation. 
For most samples the correlation function fits very well 
the exponential decay, G(r) cx e~ r l Tc . 

VI. EXPERIMENTAL RESULTS 

The parameters of studied YBCO films are presented 
in Tab. 1. The fourth column shows the width, ATr, of 
resistive transition defined as the doubled dispersion of 



the Gaussian fitting dR/dT peak for H = 0. The value 
ATr equals to approximately 0.8 of the transition width 
defined by 10%-9Q% level of normal resistance. 

The width AT C of the experimentally determined dis- 
tribution function f(T c ) was calculated by the same pro- 
cedure as ATr. For samples marked by (*) the distribu- 
tion function had two rather than one peak. In this case 
we calculated AT C as a mean-squared deviation: 

AT C = 2^((T C -T C ) 2 ) (19) 

where the averaging is performed over the area under 
the double-peak Gaussian fitting f(T c ). Application of 
Eq. ( |l9| ) to the distribution function itself is less reliable 
because the value of AT C is strongly affected by the tails 
of the distribution. 

Further, we assume that the total broadening of the 
transition is caused by summation of homogeneous and 
inhomogeneous broadening and the simple relation can 
be written: 

ATj. = AT h 2 om + AT C 2 , (20) 

where AThom is the homogeneous broadening of the tran- 
sition. 

The scale r c of T c -inhomogeneity was determined by 
fitting the correlation function G(r), Eq. (|l8|), with an 
exponential decay, exp(— r/r c ). Values of r c vary much 
for different samples and depend primarily on the sub- 
strate. This is consistent with results of x-ray studies 
which revealed clusters of dislocations of ^80 /j,m size in 
MgO substrate used for sample 1. By contrast, sample 4 
grown on NdGa03 substrate was of higher quality and 
no large-scale clusters in the substrate were observed. It 
should be noted, that values of r c in Table 1 can over- 
estimate the true correlation length of T c -inhomogeneity 
especially for samples with small r c . The reason is that 
r c is always larger than the resolution of the experimental 
method, r cxp = 2 fim. Presence of Tc-inhomogeneity on 
small scales can be revealed by x-ray diffraction studies. 



Table 1. Some characteristics of studied YBCO thin film samples. The transition width, ATr, is defined by the width 
of dR/dT peak; AT C is the dispersion of T c - distribution; ATh om is the intrinsic broadening of the transition; ATebiv 
is the average width of the local temperature dependence of EBIV; r c is the correlation length of T c -distribution; 
p n (T) is the linear fit for the temperature dependence of resistivity in the 150-300 K range. 



No. 


Substrate 


T c , K 


ATr, K 


AT C , K 


AThom — 

VAT 2 - AT 2 


ATebiv, K 


r c , jitm 


p n (T, K), /iQm 


1 


MgO 


86 


1.5 


1.2 


1.1 


1.1 


80 


1.06+0.0035 T 


2 


AlLa0 3 


92.8 


0.4 


0.3 


0.3 


0.3 




0.14+0.003 T 


3 


AlLa0 3 


91.5 


1.7 


0.8* 


1.5 


0.9 


45 




4 


NdGa0 3 


89 


1.5 


0.2 


1.5 


0.9 


6 


1.21+ 0.02 T 


5 


NdGa0 3 


88.5 


1.9 


0.3 


1.9 


1.7 


16 


1.26+0.018 T 


6 


NdGa0 3 


87 


1.7 


0.8* 


1.5 


0.5 


33 


0.96+0.0045 T 
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The size of the area where the coherent scattering of x- 
ray wave is established has been found to be 30-10(M 
for YBCO filmsJaH This value defines the lower limit 
for r c . It is in agreement with the value r c « 30 A for 
the size of T c -uniform fragment in YBCO film deduced 
from analysis of experimental data on-voltage noise in 
the superconducting transition region.cil 

The seventh column, ATebiVj shows the average width 
of the local temperature dependence of the EBIV which 
should be closely related to the homogeneous broaden- 
ing AXhom- Indeed, a good agreement is observed for 
samples 1, 2 and 5. In the case of samples 3 and 6 a 
specific shape of distribution function f(T c ) which in- 
validates simple relation ( pcf ) can be responsible for the 
deviation. For sample 4 this deviation is probably related 
to very short correlation length r c . The last column in 
Tab. 1 represents the linear fit for the temperature de- 
pendence of resistivity in the 150-300 K range; the error 
in determination of the fit coefficients is 0.2-1%. 

As follows from Tab. 1, inhomogencous broadening, 
AT C , of the resistive transition is of the same order as 
homogeneous one, AThom- The homogeneous broaden- 
ing characterizes the transition width for a fragment oi 
superconducting film of 2 /im size. This width can be 
either an intrinsic property of a homogeneous supercon- 
ductor or it can be associated with T c -inhomogeneity on 
scales < 2 /im. Large scatter of ATh om in Tab. 1 suggests 
the presence of small-scale Tc-inhomogeneity at least in 
the samples with large AX^om. 

Let us now examine the effect of Tc-inhomogeneity on 
the experimental temperature dependences of conductiv- 
ity. Data for samples 1 and 4 with maximal and minimal 
r c will be analyzed. In order to extract the supercon- 
ducting contribution, a s (T,H), to conductivity from the 
measured resistance R(T, H) we use Eq. ([l7|) and data 
from Tab. 1. The extracted temperature dependences 
of cj s were fitted by two models: for homogeneous and 
for Tc-inhomogeneous superconductor. For homogeneous 
superconductor they were fitted directly by low-field ap- 
proximation, Eqs. (^|) and (^), derived in Sec. II. The pa- 
rameters Co j K arL< i T c were free. For T c -inhomogeneous 
superconductor the same formulas were used to calculate 
conductivities of local fragments with uniform T c . Effec- 
tive conductivity a s (T, H ) of the whole sample was calcu- 
lated by solving resistor networks based on the measured 
T c -maps. This method has only two fitting parameters: 
<To and k. Theoretical estimate for ctq was obtained in 
Sec. II by comparison of the results of UD model at high 
temperatures and Aslamazov-Larkin formula. However, 
because of sample imperfections and a large error in de- 
termination of the sample thickness, (To should be a free 
parameter. For studied samples a® differs from the value 
given by Eq. ( |l3| ) by a factor between 0.6 to 1.5. As fol- 
lows from formulas of Sec. II, <7o controls the magnitude 
of the Cooper pair conductivity, while k determines the 
width of the resistive transition. 

The experimental dependences cr s (T) for sample 1 for 
three magnetic fields and their fits by the "inhomoge- 



neous" model are presented in Fig. g. The dashed line 
shows the fit by the "homogeneous" model for H = 0.6 T. 
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FIG. 6. Experimental temperature dependences (sym- 
bols) of superconducting contribution to conductivity, a s , for 
sample 1 and their fits (solid lines) calculated by solving resis- 
tor networks based on the measured spatial distribution of T c 
and Eqs. (|) and @for three magnetic fields: (1) H = 0.3 T, 
(2) H = 0.6 T, (3) H = 0.9 T; the fitting parameter is 
K = 50. The dashed line shows fit for a homogeneous su- 
perconductor, Eqs. (^) and (^|), for H = 0.6 T with fitting 
parameters T c — 87.1 K and k — 30. Fits taking account of 
T c -inhomogeneity are in a much better agreement with exper- 
iment. 
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FIG. 7. Experimental temperature dependences (sym- 
bols) of superconducting contribution to conductivity, a s , for 
sample 4 and their fits (lines) based on Eqs. (^) and (^|)with 
account of T c -inhomogeneity for three magnetic fields: (1) 
H = 0.3 T, (2) H = 0.6 T, (3) H = 0.9 T. Since the cor- 
relation length r c of T c -inhomogeneity for sample 4 is very 
small and comparable to the resolution of T c -map, the effec- 
tive medium approximation with Gaussian T c -distribution is 
used for calculations. The fitting parameters are T c = 90 K, 
and k — 30. 
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It can be seen that this model strongly deviates from 
the experimental curve. The "homogeneous" model pre- 
dicts an abrupt rise in conductivity as the temperature 
decreases which is not observed on experiment. In "in- 
homogeneous" model this contradiction disappears. 

As it can be seen from Tab. 1, the width, AT C , of 
the measured T c -distribution for sample 4 is substan- 
tially less than the transition width, ATr. We believe 
that this fact as well as small r c are related to presence 
of T c -inhomogeneities on scales less than the experimen- 
tal resolution, r oxp . In this case calculations based on the 
measured T c -map are not reliable. Instead, in order to 
calculate the effective conductivity for sample 4, we used 
EMA and a Gaussian T c -distribution function. The re- 
sults for three magnetic fields are shown in Fig. Addi- 
tional fitting parameters, the average and the dispersion 
of Gaussian distribution, were found to be T c q = 89.1 K, 
and 6T C = 1.3 K. 

Presence of small-scale T c -inhomogeneities is probably 
the reason for difference in n determined from fitting the 
experimental u s (T) dependences for different samples. 
For sample 4 (Fig. 0) the best fits based on the EMA 
are obtained with k = 30, while for sample 1 (Fig. ||) the 
best fits based on the measured T c -map give k = 50. The 
high k in the latter case leads to additional broadening 
of the transition compensating lack of information about 
small-scale T c -inhomogeneities. Thus, the value k = 30 
is more reliable and it is used for calculations presented 
in Sec. III. 

To summarize, there arc two ways to take T c - 
inhomogeneity into account: direct resistor network cal- 
culations based on T c -map, and EMA along with a Gaus- 
sian T c -distribution. The resistor network calculations 
have the advantage of using actual spatial distribution 
of T c in the sample. It has the information about loca- 
tion of regions with various T c allowing the calculation of 
percolative current distribution in given HTSC film. On 
the other hand, the drawback of this model is that the 
T c -map is measured with finite spatial resolution. Thus, 
one should use either T c -map or EMA for large and small 
values of the correlation length, r c , respectively. 

In Fig. [t] all models significantly deviate from the ex- 
perimental data at sufficiently low temperatures. We ex- 
plain this deviation by the vortex pinning which comes 
into play for low temperatures and prevents the dissipa- 
tion associated with flux flow. UD model does not take 
the pinning into account and, thus, overestimates the 
dissipation rate. We believe that in the low-temperature 
part of the superconducting transition it is the strength 
and concentration of pinning centers rather than Tc- 
distribution that controls the transport properties. 

It is well-known that while the resistance of an inho- 
mogeneous system is determined by the second moment 
of current distribution, the resistance fluctuations are de- 
termined by the fourth moment. Therefore the resistance 
noise is far more sensitive to the presence-of all kind of 
inhomogeneities than the resistance itself.Ej This means 
that although this work presents analysis of the trans- 



port properties only, one can expect far stronger effect of 
T c -inhomogeneity on the noise properties of superconduc- 
tors. Even simple analysis not involving any particular 
dependence of local conductivity on T and H shows a 
strong effect oLT c -inhomogeneity on the level of thermo- 
dynamic noiseEiand the noise associated with fluctua- 
tions of local T c £3 

The properties of single crystals differ much from those 
of thin films and need a special consideration. It is gen- 
erally believed that small transition width, 0.1 — 0.3 K, 
in zero magnetic field observed in single crystals proves 
their high homogeneity. Therefore, the experimental 
data on single crystals are often used to get insight 
into fundamental intrinsic properties of superconduc- 
tors. Nevertheless, recent theoretical and experimental 
investigations make their homogeneity, in garticular, T c - 
homogeneity questionable. It is predictedQ that various 
extended structural defects, e.g., dislocations can give 
rise to formation of the extended regions with enhanced 
T c nearby. Studies of the influence of oxygen stoichiom- 
etry on the magnetization curves of YBCO crystals sug- 
gest that the so called peak effect widely observed in 
HTSC crystals is associated with the presence of local 
regions with— reduced oxygen content, and, hence, re- 
duced T c .Ei]0 Presence of non-uniform T c -distribution 
in YBa2Cu307_5 and Bi 2 Sr 2 CaCu208 crystals follows 
from experimental data on in-plane magnetoresistivity 
anomalitiesJl3 Further, large-scale spatial variations of 
oxygen composition, implying variations of T c , were 
observedcl in YBCO single crystals by x-ray studies. 
However, the spatial scale of Tc-inhomogeneities in crys- 
tals often- has a value comparable to the size of the 
sample BES In such a case, despite a wide distribution of 
T c over the sample, the superconducting transition can 
be very sharp because of a percolation over high-T c re- 
gions along one of the sample edges. Unfortunately, such 
situations cannot be properly treated in the frame of the 
effective medium approach because it assumes purely un- 
corrected T;-distribution. EMA can neither be applica- 
ble to describe wires of higher T c near extended structural 
defectsE Thus, we do not expect that the results of this 
work would be applicable to HTSC crystals. Neverthe- 
less, there are grounds to believe that inhomogeneity of 
crystals strongly manifests itself in their properties and 
deserves a detailed analysis. 



VII. CONCLUSIONS 

T c -inhomogeneity of YBCO films is directly demon- 
strated by measuring spatial distributions of T c by low- 
temperature SEM with 2 /im resolution. The dispersion 
of T;-distribution was found to be of the order of 1 K 
which is comparable to the resistive transition width. 
This result indicates inhomogeneous broadening of the 
resistive transition for the films studied. 

We obtain a non-explicit expression for Cooper pair 
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conductivity a s (T, H,T C ) of a homogeneous supercon- 
ductor, which is valid throughout the transition region 
for magnetic fields H <C H c2 {Q). For YBa 2 Cu307_(5, 
it can be reduced to an explicit expression for fields 

h «o.m c2 (o). 

We find that the error in the apparent value of 
cr s (T, H, T c ) due to T c -inhomogeneity is maximal for low 
magnetic fields and temperatures close to T c . For YBCO 
films with a Gaussian T c -distribution with 1 K dispersion, 
ignoring T c -inhomogeneity leads to more than 30% error 
in a s in the region restricted to temperatures \T — T c \ < 
0.5 K and magnetic fields H < 1 T. Thus, it is necessary 
to be cautious when carrying out quantitative analysis 
of experimental data in the transition region. One of the 
following is recommended: (i) carry out all measurements 
beyond the region affected by T c -inhomogeneity, i.e., at 
very high magnetic fields or at temperatures far from T c ; 
(ii) take T c -inhomogeneity into account by measuring T c - 
spatial distribution or, at least, by assuming a Gaussian 
distribution and using EMA or similar approximation. 

Finally, it should be noted that the boundaries of H- 
T plane region affected by T c -inhomogeneity are deter- 
mined not only by microscopic superconducting parame- 
ters, but also by material parameters such as dispersion 
and correlation length of T c -inhomogeneity. Neverthe- 
less, a transition width of the order of 1 K seems typical 
for YBCO films, while Bi-bascd films usually have even 
broader transition. Thus, the presented results are likely 
to be relevant to most HTSC films. 
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